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ABSTRACT 

The magnetorotational instability (MRI) of a weakly ionized, differentially rotat- 
ing, magnetized plasma disc is investigated in the multi-fluid framework. The disc 
is threaded by a uniform vertical magnetic field and charge is carried by electrons and 
ions only. The inclusion of ion inertia causes significant modification to the conduc- 
tivity tensor in a weakly ionized disc. The parallel, Pedersen and Hall component of 
conductivity tensor becomes time dependent quantities resulting in ac and dc compo- 
nents of the conductivity. The time dependence of the conductivity causes significant 
modification to the parameter window of magnetorotational instability . 

The effect of ambipolar and Hall diffusion on the linear growth of the magnetoro- 
tational instability is examined in the presence of time dependent conductivity tensor. 
We find that the growth rate in the ambipolar regime can become somewhat larger 
than the rotational frequency, especially when the departure from ideal MHD is sig- 
nificant. Further, the instability operates on large scale lengths. This has important 
implication for angular momentum transport in the disc. 

When charged grains are the dominant ions, their inertia will play important role 
near the mid plane of the protoplanetary discs. Ion inertia could also be important 
in transporting angular momentum in accretion discs around compact objects, in cat- 
aclysmic variables. For example, in cataclysmic variables, where mass flows from a 
companion main sequence star on to a white dwarf, the ionization fraction in the disc 
can vary in a wide range. The ion inertial effect in such a disc could significantly mod- 
ify the magnetorotational instability and threfore, this instability could be a possible 
driver of the observed turbulent motion. 

Key words: magnetohydrodynamics, star formation, accretion discs, charged grains, 
magnetorotational instability . 



1 INTRODUCTION 



Angular momentuni transport has long been recognised as 
a key issue in accretion disc theories (Lynden-Bell 1969; 
Sakura & Sunyaev 1973). However, until the 1990s, a viable 
physical mechanism necessary to facilitate this transport in 
the absence of tidal effects or gravitational instabilities was 
unknown. The Balbus-Hawley (or magnetorotational) insta- 
bility (Velikhov 1959; Chandrasekhar 1961) was proposed 
(Balbus & Hawley 1991; Hawley&Balbus 1992; Balbus & 
Hawley 1998) as a viable mechanism that can efficiently 
drive MHD turbulence and transport angular momentum in 
the disc. This opened the door for its application to a wide 
variety of astrophysical discs. The requirement for the mag- 
netorotational instability to operate in such a disc is that 
the ambient magnetic field is subthermal at the disc mid- 
plane and is well coupled to the disc matter. Although, the 
lower bound on the weak, subthermal field has never been 



specified, in recent work this issue has been addressed in 
the framework of fully ionized, coUisionless cold electron-ion 
plasma (Krolik & Zweibel 2006). For a highly ionized disc, 
the requirement of a weak, subthermal field is easily satis- 
fied and the magnetorotational instability grows at the rota- 
tion frequency Q. of the disc as a low frequency Alfven mode 
with fc Va ~ fi, where k is the wavenumber and Va is the 
Alfven velocity. However, many astrophysical discs are not 
well coupled to the magnetic field. Circumstellar, Protoplan- 
etary (PPD), Dwarf Novae (DN), and, proto-neutron-star 
discs are good examples of weakly ionized discs with very 
low to low (PPDs) and high (DNs) fractional ionization. In 
PPDs for example, the sources of ionization are limited to 
the disc surface and magnetorotational instability may op- 
erate only in the outer envelope of the disc (Gammie 1996) 
unless some nonthermal source of ionization viz. the colli- 
sion of the energetic electrons with neutrals or magnetorota- 
tional instability induced turbulent convective homogeniza- 
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tion of the entire disc (Inutsuka & Sano 2005) is assumed. 
DN discs are thought to have both hot and fully ionized 
accretion state as well as cold and mostly neutral accretion 
states (Cannizzo 1993; Gammie & Menou 1998). Therefore, 
the direct application of Balbus & Hawley (1991) results are 
difhcult in a weakly ionized disc. 

The effect of non-ideal MHD on the magnetorotational 
instability has been investigated by several authors: in the 
ambipolar regime (Blaes & Balbus 1994), hereafter BB94, 
(Maclow et al. 1995; Hawley & Stone 1998; Kunz & Balbus 
2004), the resistive regime (Jin 1996; Papaloizou & Terquem 
1997; Balbus & Hawley 1998; Sano et al. 1998; Sano & 
Miyama 1999; Fleming et al. 2000; Sano et al. 2000; Stone 
& Fleming 2003) and the Hall regime (Wardle 1999) W99 
hereafter; (Balbus & Terquem 2001), BTOl hereafter; (Sano 
& Stone 2002a,b; Salmeron & Wardle 2003, 2005; Desch 
2004). At the densities relevant to cloud cores, ambipolar 
and Hall diffusion plays an important role in the transport 
of mass and angular momentum (Wardle & Ng 1999; Bal- 
bus & Terquem 2001). W99 and BTOl found that collision 
of neutrals with the ionized gas in a weakly ionized disc 
determines the relative importance of Ambipolar, Hall or 
Ohmic diffusion on the magnetorotational instability . The 
ambipolar and Hall effects are particularly important when 
the ionization in the disc is very low and the departure from 
ideal MHD is severe. 

The dynamics of a weakly ionized disc was investigated 
in the limit of zero inertia of the ionized plasma components 
by W99 and BTOl. This is usually an excellent approxima- 
tion when the fractional ionization is low, and allows the ion- 
ized components of the fluid (viz electrons, ions and grains) 
to be treated on an equal footing. However, there are situ- 
ations - even in the low fractional ionisation limit - where 
the inertia of the charged species is important in determining 
their drift with respect to the neutral component and hence 
the diffusion of the magnetic field. In the weakly-ionised 
limit this becomes important when the inertial terms in the 
ion equation of motion start to compete with the magnetic 
and neutral collision terms, in other words when the disc 
frequency Q. becomes comparable to or exceeds the collision 
frequency with neutrals and the gyrofrequency. This effect 
may become important for charged dust grains because their 
large mass implies low collision and gyrofrequencies. For ex- 
ample, in dense PPDs, the dominant ion species are posi- 
tively charged grains (especially when ~ uh ^ 10^^ cm^'*), 
(Wardle & Ng 1999). In the high fractional ionization re- 
gions also, e.g. near the surface of a magnetic cataclysmic 
variable, ion inertia may become important (Warner 1995). 

The effect of ion-inertia on the magnetorotational in- 
stability in a two-fiuid framework was considered by BB94. 
The magnetic flux was assumed frozen into the plasma com- 
ponent in their formulation. This limits the applicabilty of 
their results to the ambipolar diffusion regime. We know 
from W99 and BTOl that Hall effect can compete with am- 
bipolar diffusion in the weakly ionized regions of the disc. In 
the present work we adopt a three component model - neu- 
trals, ions and electrons where by "ion" and "electron" we 
mean the most massive and least massive charged species, 
whether positively or negatively charged - to show that ion 
inertial effects modify the growth rate substantially and in- 
creases the parameter window in which the instability may 
operate. In particular, the growth rate in the presence of ion 



inertia may exceed the Oort A value (0.75 fi) limit. Growth 
rates larger than the Oort A limit have recently been re- 
ported for collisionless plasmas in the presence of plasma 
kinetic/viscous effects (Quataert et al. 2002; Sharma et al. 
2003; Balbus 2004), where kinetic and MHD effects com- 
bine to give a high growth rate (1.7 il) and shift the fastest 
growing modes towards longer wavelengths. 

We shall give a general formulation of the problem that 
will not only cover the regions of applicability of BB94 and 
W99 but also cover the unexplored regions. This paper revis- 
its the magnetorotational instability in a weakly ionized disc 
(W99) by incorporating the effect of ion inertia. In section 2, 
we discuss formulation of BB94 and compare and contrast 
the region of applicability of BB94, W99 in the context of 
present work. In section 3, a general formulation of a weakly 
ionized, near-Keplerian, magnetized disc is given. In section 
4 we describe equilibrium and linearization of the disc and 
derive the dispersion relation. Section 5 discusses the ener- 
getics of a weakly ionized, magnetized disc. In this section, 
we first discuss the role played by the ambipolar and Hall 
terms in the neutral dynamics. Further, using energy ar- 
guments, we discuss the conditions under which ambipolar 
diffusion can proceed without dissipation and when ambipo- 
lar diffusion can destabilize the disc. We also discuss how 
theHall term can destabilize the magnetorotational instabil- 
ity . The properties of wave helicity, is derived. In section 
6, we give the detailed numerical solution of the dispersion 
relation. In section 7, application of the result to various as- 
trophysical discs is discussed. Section 8 presents a summary 
of our results. 



2 ION INERTIA AND BB94 

BB94 adopted the ambipolar diffusion approximation, in 
which the magnetic field is frozen into an electrically neutral 
ionised plasma coupled by collisions to the neutral fluid. 

The effect of ion inertia on the dynamics [pi dtVi where 
pi is the ion mass density and dt = d/dt = dt + Vi ■ V 
is the convective derivative), is retained in the momentum 
equation for the ionised component but its effect on Ohm's 
law is ignored in BB94. This limits the applicability of the 
results to the ambipolar diffusion limit. 

To better appreciate this point, let us briefly recast the 
two fluid formulation of BB94 starting with seperate ion 
and electron fluid equations. The equation of motion for the 
ionised fluid is derived assuming that magnetic fleld is frozen 
in the electron fluid. 



= —ene{E + Ve X B/c) 



(1) 



and summing the electron and ion momentum equations, to 
yield 

dv^ VPi , _ , , . . J X B _, 

dt Pi c Pi 

equation (15) of BB94. Here e is the electronic charge, n^ is 
the electron number density, Vin is the ion-neutral collision 
frequency, Ve, Vi and v„ are the electron, ion and neutral 
velocities, $ is the gravitational potential, Pi is the ion pres- 
sure, E, B are the electric and magnetic flelds and c is the 
speed of light and J — ene{vi — Ve) is the current den- 
sity. Taking the curl of (1) and using Maxwell's equation. 
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c\/xE = -dtB, will give dtB ^ \/ X {v^ X B), i.e. the 
magnetic field is convected away by the electron fluid. If we 
want to express the right hand side of induction equation in 
terms of ion velocity, we obtain 



^=V X (u, X B)-VxFh, 

where the Hall term is 
J xB 



Fh^ 



(4) 



Since J X B/c ~ pi {dtVi + Ui„Vi) (here Ui„ is 
the ion-neutral collision frequency), the Hall term can be 
dropped from the induction equation only if {vin,uj) <§C 
i^ci {= e B / {mi c)) , i.e. the ion-gyration period (io^^) is 
smaller/ faster than the dynamical time (iu^^) or the ion- 
neutral collision time {v'^ ). We see that unlike a two com- 
ponent electron-ion plasma, where the Hall term can be in- 
troduced only through the ion inertial term (dtVi), in a 
weakly ionized multi-component plasma, the Hall effect ap- 
pears either via ion-neutral collision or via the dt Vi term or 
both. 

Replacing dt B by AB/St and VxFh by 
cB AB/{4:Trnieh Ax), one sees that the Hall term 



scales as l/{ujpiL), (here tUp 



(4 7rc rii/mij ' is the 



ion plasma frequency and L is the characteristic size of 
the system), i.e. the Hall term is important on a scale 
shorter than the ion-inertial scale. Clearly, Hall MHD 
introduces two disparate, interacting scales, a microscopic 
scale, i.e. the ion-skin depth {Si = c/copi = VA/i-Od) and 
a macroscopic scale, the disc size. The Hall term can be 
dropped if ion-inertial effects are unimportant. Leaving out 
the effect of inertia in the induction equation but retaining 
them in dynamics is not consistent and, in such a scenario, 
one would expect that magnetorotational instability will 
merely shift towards long wavelength, as has already been 
noted by BB94. 

If we start with the ion equation of motion (BB94 (15)), 
and express the electric field as 



E 



ViX B 



■mi 
e 



dv, ^ ^^ VPi ^ 

— - -I- V* -I h Vin Vd 

dt pi 



(5) 



where vd — Vi — w„, then taking the curl and using 
Maxwell's equation, one arrives at the following induction 
equation 



dB 

~dt 



= \7x{vi X B) 



rUi c 
e 



VinVXVD +VX 



dvi 
dt 



(6) 



Here uniform density is assumed while operating with the 
curl on equation (5). This equation has an additional term 
in comparision with equation (16) of BB94, with important 
consequences on the magnetic diffusivity, since the rate of 
change of magnetic flux is given as 

d_ 
di ., ^ 

'dB 



— Bds = 



dB , 



Vi X dl- B 



Vx {v, X B) 



ds. 



(7) 



Making use of equation (6) in (7) , we get 



^ B 



ds — 

ft 



Vd 



— VX 



9vi ._ ^ 

-TT- -V^X {\/XV^) 



dl 



ds 



(8) 



Here B = B/B and use has been made of dvi/dt — dtVi 



Vi X {\7xvi)+\/vi/2. The ion HaU parameter l3i 



if l^in 



gives the ratio between the ion-cyclotron to ion-neutral col- 
lision frequencies. The above equation can be rewritten as 



(3) ^ 



1 „ 
— Vxvi 



''^ = ^i 



Vd 



dl. 



(9) 



We see from equation (9) that the generalized flux that is 
a combination of magnetic flux and vorticity is not con- 
served. The rate at which this flux decays is directly related 
to the coUisional coupling between ions and neutrals. If the 
ion magnetization is weak, i.e. the ion-cyclotron frequency is 
less than ion-neutral collision frequency {f3i -^ 0), then the 
flux-decay rate could be very large for a finite ion-neutral 
drift speed vd- However, if the relative ion- neutral drift is 
negligible, the generalized fiux is conserved irrespective of 
the ion magnetization level. Thus it is the combination of 
the magnetic flux and the vorticity that is conserved in the 
absence of collision {j3i — > cx)). The BB94 formulation as- 
sumes that magnetic flux is frozen in the ion-fluid which is 
valid if apart from ignoring the right hand side of equation 
(9), we also assume that ujci ^ oo. In this limit however, the 
role of ion inertia becomes increasingly unimportant and we 
approach W99 limit. Clearly, BB94 does not treat the ef- 
fect of ion inertia in a consistent fashion and their results 
are not applicable in most of the weakly ionized parameter 
space where the ion Hall parameter is ~ 1. In Fig.l, we plot 
the range of applicability of BB94 and W99. We see that 
BB94 is applicable when Pi ^ 1 for ujci > fi for arbitrary 
relation between Vin and Q. BTOl show that ft ^ 1 implies 
that the Hall term dominates over ambipolar term. W99 is 
applicable for uJci > ^ and Vin > Q for arbitrary l3i. There- 
fore, uJci/0. < 1 and Vin/O. < 1 is an unexplored region in 
BB94 and W99 framework. We see that in a PPD, for a mil- 
liGauss field, for a positive grains of mass 10~^^ — 10~^^ g, 
uJci/^ < 1. As has been noted elsewhere, near the mid- 
plane of PPDs, dust grains can be the dominant charged 
constituent over extended regions (~ 1 — 5 AU). For sub- 
micron sized grains 0.1 ^m, negatively charged grains domi- 
nate whenever n„ > 10^^ cm""^ and positively charged grains 
dominate for n„ > lO^"* cm"^ (Wardle & Ng 1999). Depend- 
ing upon neutral density in the disc, the ratio Vin/Q. can have 
any value and thus, it is importnat to extend the BB94 anal- 
ysis to the unexplored regions with full. Hall and ambipolar 
effects in the spirit of W99. We shall give a general deriva- 
tion of the induction equation in a weakly ionized medium 
in the text though, we adopt the conductivity approach of 
W99 to investigate the effect of ion inertia on the magne- 
torotational instability and find a significant increase in the 
growth rate. Further, the parameter window in which insta- 
bility operates expands considerably to large scale lengths. 
It will have an important bearing on the angular momentum 
transport and onset of turbulence. 



3 FORMULATION - MHD EQUATIONS 

The dynamics of a weakly ionized disc, consisting of elec- 
trons, ions, neutrals and charged and neutral dust grains, in 
the presence of a gravitational field $ of a central mass point 
M is described by a set of multifiuid equations. A multi-fluid 
approach representing each and every particle species is not 
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Figure 1. The region of applicability of BB94 and W99. 



very fruitful since depending upon the fractional ionization, 
the presence of some of the ionized component in the disc 
can be neglected. We shall assume that the ion density in 
the disc is mainly due to the presence of positively charged 
grains. Such a situation will correspond to a very dense re- 
gion of PPDs (Wardle & Ng 1999). 

As the disc matter is weakly ionized, generally the 
inertia of the charged species are neglected (W99). How- 
ever, apart from pure scientific curiosity about the effect of 
charged particle inertia on the magnetorotational instabil- 
ity , there are astrophysical environments where the plasma 
inertial effects may become important. For example, ion in- 
ertial effect may compete with coUisional and electromag- 
netic effects in the dense region of PPDs as well as in the 
high fractional ionization discs around magnetic cataclysmic 
variables. An estimate of the ion inertial, coUisional and 
electromagnetic effects are given below. Since we wish to 
investigate the role of ion dynamics on the magnetorota- 
tional instability , we shall retain the inertia term in the ion 
equation of motion. This will result in conductivity tensor 
a becoming frequency dependent. 

The basic set of equations describing a partially ionized, 
non-self-gravitating magnetized plasma disc, consisting of 
electrons, ions, and neutrals, are as follows. The continuity 
equation is 

^+V.{p,v,) = Q. (10) 

Here pj is the mass density and Vj is the velocity of the 
various plasma components. 

The momentum equations for electrons, ions and neu- 
trals are 



0: 



E' + 



Ve X B 



pe l^en "^e j 



(11) 



pidtVi = eui [ E' ^ ^- 1 — Pi Uin Vi — pidt v„, (12) 



(13) 



Pn dt Vn = - VP — p„ V<1> + ^ pj Vjn Vj 



Here Vj is the drift velocity through the neutrals, Fj = 
qj rij (E' + Vj X B/c) is the Lorentz force, rij is the num- 



ber density, and, j stands for electrons (q^ = — e), and ions 
(qi = e). Grains are assumed to have single positive electic 
charge. The electric field E' = E + Vn X B is written in 
the frame comoving with the neutrals. In equation (13), the 
gravitational potential $ due to central mass, is given by 



$ = -- 



GM 



(14) 



At the disc midplane, the Keplerian centripetal accelera- 
tion v^/r = GM/r^ balances the radial component of the 
gravitational potential. Equations (11)- (12) does not have a 
pressure gradient term since pressure effects will be negli- 
gible in a weakly ionized disc. The effects of ionization and 
recombination are also omitted from the neutral dynamics 
for the same reason. 

The collision frequency is 

^iri = "Ijn pn = \ pn- (15) 

mn+nij 
Here < av >j is the rate coefRcient for the momentum trans- 
fer by the collision of the j particle with the neutrals. The 
ion-neutral and electron-neutral rate coefficients are (Draine 
et al. 1983) 



< cr-u >i„= 1.9 X 10 ' cm^s ^ 



< av >en= 4.5 X 10 



f T 



V 30K 



(16) 



Adopting a value of nii — 30 nip for ion mass and m„ = 
2.33 rrip for mean neutral mass where rxip — 1.67 x 10"^"* g 
is the proton mass, the ion neutral collision frequency can 
be written as 

Vrn=Pnl= ; = 1.4X10 n„ S . (17) 

rrii + m„ 

This also gives the limiting value for very small grains 
(~ 3 — 3000 ry4). For larger, micron sized grains, ion-neutral 
collision rate can vary between 10"^" to 10~^ for sizes rang- 
ing between a few Angstrom to a few microns. This can be 
seen if we write the collision rate as (Nakano & Umebayashi 
1986) 



< aV >in- 



(18) 



where T30 is the gas temperature and a_5 is the grain radius 
in units of 30 K and 10~^ cm respectively. 

We shall rewrite equations (10)-(13) in the local Keple- 
rian frame. Thus, velocity v represents the departure from 
the Keplerian motion; the fluid velocity in the laboratory 
frame is i; -f vk and dt is dt + Q 84, in the laboratory frame, 
where vk = \/G M/r (f> is the Keplerian velocity in the 
canonical cylindrical coordinate system {r,(j),z). 

Noting that near the disc midplane, on scales small com- 
pared with the disc thickness, the radial gradient in gravita- 
tional potential will be exactly cancelled by the centripetal 
term due to Keplerian motion, and, (r, <j)) component of the 
equation (12)-(13) in the absence of any tidal effects, can be 
rewritten as, 



Avi 



e 
rUi 



E'- 



AVn = - 



VPn 



Vi X B 



J X B 

CPn 



■,Vr - Av„, 



.Pn 



(19) 



(20) 
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dt 

0.5 n 



-2Q 

dt 



where operator A = 

equation can be written as 

dtB = Vx (t; X B)-cVxE' -1.5QB. 



The induction 



(21) 



In (21), VxE' contains the effect of non-ideal MHD and 
the last term accounts for the generation of the toroidal 
field from the poloidal one due to differential rotation of the 
disc (W99). 

Before describing the conductivity approach and giv- 
ing a formulation to the problem at hand, we shall give a 
general derivation of the induction equation starting from 
equation (11) and estimate the range of applicability of 
the ion-inertial effects in the accretion discs. Writing E' — 
— Ue X B/c + rj J as 



E' = 

Here 

c 
77 = 



Vi X B 



+ riJ + Fh. 



rrie n„ < a v >e 



An Tie e^ 



4-K 



(22) 



(23) 



is the electrical resistivity of the gas and Fh — J X B/eUe 
is the Hall term. Even in the absence of ion inertia - the 
so called strong coupling approximation, (Shu 1983), the 
ambipolar term modifies the induction equation, i.e. when 
piVinVi = J X B/c, 



E'^f)J + FH- 



(J X B)x B 



C pi Vin 

When ion inertial terms are also present. 



J X B 



Vt = 



-LiHi + o ( — — 

Vin dt \Pn' pn 



C pi Vin 

and, the generalized Ohm's law becomes 
(J X B/c) 



E' = riJ + FH 



Pi Vi, 



- AiVi 



X B, 



(24) 



(25) 



(26) 



where Ai — v^^ A. Now taking the curl of equation (26), 
one may write the generalized induction equation as 



dB 

~dt 



Vx 



V X B 



A-K-qJ J X B 



+Vx 



{J X B)xB 



C pi Vin 



(AiVi) X B 



(27) 



The induction equation (27) on the right hand side contains 
inductive, Ohmic diffusion. Hall, ambipolar, and ion inertial 
terms respectively. The set of equation (10), (12), (13) and 
(27) can be closed by an equation of state. 

Assuming that ion-inertial time scale is of the order of 
disc rotation period, we may write Ai ~ Q./vin. Then the 
ratio of the ion inertial term to the inductive term can be 
written as 



\v X B\ 
\A-i_Vi X B\ 






(28) 



A different ratio that measures the coupling of the neutral 
to the ion with respect to the Keplerian frequency have been 
used by BB94 (Menou & Quataert 2001) 



ly-fii i^tn 



n 



fi 



(29) 



where a — pi/ pn. It is clear from equation (29) that magne- 
torotational instability can act on both the neutral as well 
as on the ion fluid simultaneously if a ~ 1. 

Assuming an equation of state or dropping the pressure 
gradient term in the neutral equation of motion, equations 
(10), (19), (20) and (27) with Maxwell's equations 



V X B = — J, 

c 



V • B = 0, 



(30) 



form a complete set. 



4 LINEARIZATION 

We consider a thin disc implying that the radial scale over 
which physical quantities vary is much larger than the disc 
scale height, H — C's/^. The initial steady state is assumed 
uniform and homogeneous with a vertical magnetic field 
B — Bz and zero v, w„, VP, E' and J. 

We shall assume transverse fluctuations and denote re- 
sulting two-dimensional vectors by subscript _L, to investi- 
gate Alfven modes in the disc. We seek plane wave solution 
of the form exp i{ujt ~ kz). 



4.1 The conductivity tensor 

The conductivity tensor (t can be found by considering the 
drifts of charged particles in response to the electromagnetic 
field (Cowling 1957; Norman & Heyvaerts 1985; Nakano & 
Umebayashi 1986; Wardle & Ng 1999). We shaU first derive 
conductivity tensor from equation of motion for charged par- 
ticles, from (11) and (19) by eliminating Ve and Vi in favour 
of E' and B and then give the expressions for parallel. Hall 
and Pedersen component of this tensor o in the generalized 
Ohm's law, J = cr E' . 

Since we assume a homogeneous steady state, it implies 
that the inohomegenity scale length (in this case verticle 
scale height H of the disk) is much larger than the charac- 
teristice wavelength of the normal modes in the disc. Thus, 
the contribution of the pressure gradient will be neglected 
while inverting equation (13). Furthermore, we shall also ig- 
nore terms of the order ~ pi/ pn in the conductivity tensor. 
Expressing velocities Vj in equations (11) and (19) in terms 
of electric field E' the relationship between J and E' can 
be written as, 



J ^(T E' = a«E'u + anB X E'j_ + apE'j_ 



(31) 



where a|| , an and ap are the field-parallel. Hall and Peder- 
sen components of the conductivity tensor a and B is the 
unit vector along the magnetic field. If charged species j has 
particle mass rrij, charge ZjC, number density Uj, then the 
Hall parameter is given as 

ZjcB m -\~ rrij 



/3.— 



rrijC < av >j p 



(32) 



where m is the mean neutral particle mass and have dropped 
subscript 'n' from neutral quantities. As noted in section 2, 
the Hall parameter determines the magnitude of the mag- 
netic flux transport. The ratio of ion to electron Hall param- 
eter suggest that in the protostellar discs, /3i//3e ~ 10~^ ^ 
1. Recall that BB94 formulation is valid when the ion-Hall 
parameter is large {f3i ^ 1), i.e. the when Lorentz force 
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dominates the ion-neutral coUisional momentum exchange. 
However, this hmit implies strongly magnetized ions and in- 
frequent collisions. 

The conductivity tensor is frequency dependent in the 
presence of ion inertia. The parallel conductivity is 



ecui 
B 



I3e + 



ft 



1 + 



-jft 



1 + 



(33) 



Since the plasma is quasi- neutral, we have assumed rie = rn. 
The conductivity has become complex. In the low frequency 
limit, when ion inertial effect is unimportant, i.e. u)/vin -^ 0, 
(T|l reduces to W99. In tu/vin -^ oo, o-|| has both a real and 
an imaginary component 



i?e[cr||] ~ —p^-Pe, 



B 



and. 



7m[cr||] 



ecUi / uj 



Gp ^ Op 



B V uj 
The Pedersen conductivity is 
Ao-p" 



1 + 



xO 



Here frequency independent part, op is 



e c Y^ Uj 



-/3 



and frequency dependent part, Aap is 



Ao-p 



B 



-/3? 



(OM-1). 



(34) 
(35) 

(36) 

(37) 

(38) 



Here Q{ij) = {1 + f3f) D1/D2, Di- ito/uin + 1, D2 = DJ + 
Oi CI2, Cli ^ 2C1 + Pi, f22 = 0.5f2 -I- Pi and Cl = Q,/vin. In 
order to isolate real and imaginary part of ap and investigate 
the low and high frequency limits, we write Re[Q{ijj)] and 
Im[Q{uj)\'ds 



(39) 



Re[Q{Lo)] 




l+( 


.^j +" 


Lf^2 




1 + /3? 


l-(l^)' + ^i"2 


\u,f 


and. 


Im[Q{Lj)] 
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i + (t^)'-"i^2 




1+/3? 


1 { " 
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^^(^ 


.)' 



(40) 



In the low frequency limit, uj/i^in -^ 0, and assuming cj ~ il, 
/?e[Q(cj)] ~ 1 and Im[Q{uj)] ~ 0. Thus ap = cr?.. In the 
high frequency limit, when uj/i^in -^ 00, Re[Q{uj)] ~ and 
Im[Q{Lo)] « —Vin/uj. Thus, 

ec Weft 



i?e[o-p] = 



and. 



B 1+/3I 



7m[ap] = — —n^Pi { 

B \ Uin 



(41) 



(42) 



Like the parallel conductivity, the real part of the Peder- 
sen conductivity in the high frequency limit is mainly due 



to electron magnetization and imaginary part is due to ion 
magnetization. Complex resistivity is well known in LCR 
circuits where the resonance condition is found by setting 
imaginary part of the impedance to zero. In the present case, 
a similar resonance condition can be derived by setting nu- 
merator of equation (40) to zero, i.e. cu^ /uf„ = —1 + Qi 0,2, 
we get 



1 + 



Ao-p 



^1^2 



2 l-f. 



(43) 



It is important to note that the scale of frequency depen- 
dent conductivity associated with the resonance {iuj ~ Q), 
can become larger than the dc conductivity. We anticipate, 
therefore, that the frequency dependent Pedersen conduc- 
tivity will significantly modify the magnetorotational insta- 
bility . 

The Hall conductivity is 



CTH ^ — 



1 


1 


„<L 



„(1 



where the dc part is given as 



ecUi 



Et^^ 



+ /3? 



(44) 



(45) 



The frequency dependent parts, Aanr, and, (Th4> are 
e c riiPf 

'BTTpf 



AaHr=^-^i^{H,H-l) 



AaHd, = 



and. 



Hjico) = 



e c riiPi 



B l+Pi 



■{H2{CU) - 1) , 



(i + p?) n, 



(46) 



(47) 



P^D2 ' 

for j = 1,2. The radial and azimuthal component of the Hall 
conductivities are not equal due to the unequal radial and 
azimuthal coefficient in the ion momentum equation. The 
real and imaginary part of Hj is given by equations (39)- 
(40) if we recognize that right hand side of equation (47) has 
a factor I/-D2 ~ Q{cj)/Di (1 + pf). Thus, near resonance 



Aoj/r 



Qi 1 1 + 010-2 



ft 



l-h^ir^a +4 



(48) 



Except for fii becoming f22, the remainder of the expression 
for the azimuthal factor will be identical to equation (48). 

The analogy to LCR resonance can be brought closer if 
we express Ohm's law, equation (31) in diagonal form. To 
that end, we shall express E'j_ in the eigen-basis vectors of 
the rotation operator e± = (e^, ±iej,)/\/2. Then B X e± — 
=pi e^ and, Ohm's law for the transverse component J± can 
be written as 



J± = {crp =F icrn) E'j_ 



(49) 



While writing equation (49) , we have assumed anr = o-h^ ~ 
an- We may write J± = Z± (tu) E'j__^. Here the impedance 
Z± (uj) — R± + i X± (lo), with R± = Re[ap] ± Im[aH], and, 
X± — Im[ap] =F Re[aH]- Near resonance, Z± = R±. Thus 
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R±E'^ 



El^ 



1 + 



For a growing mode ioj ~ f2, and, oscillations in the current 
can take place in the absence of a neutral-frame electric field 
E'j_ — + 0. In a weakly ionized disc, if /3i < 1, the oscillation 
in the current is set by the ions diffusing across the ambient 
magnetic field. For Vin ^ iui, near resonance, collision will 
act like a driver of the resonance. 

We note that near resonance, the conductivity may 
change sign. From (43), when il < Vin, the negative conduc- 
tivity will play an important role. The ratio of the dynami- 
cal to the ion-neutral collision frequency determines whether 
negative conductivity is important. The DC conductivity be- 
coming negative within certain frequency range, in the mi- 
crowave irradiation is well known in the condensed matter 
literature, e.g. Ryzhii (2005). 

4.2 Dispersion relation 

The linearized neutral equation of motion (20) can be writ- 
ten as 



-in 

2 



2iQ 



Sv± — 



I 2 
-kVA 






(51) 



where subscript _L denotes two dimensional vector in the disc 
plane and va ~ B/y^An p is the Alfven velocity in the total 
fluid. The linearized induction equation, after substituting 
for dv± is given as 



f cj'i + 3^2 


2iujn 
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ly -2iujD. 
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where uj\ = 


LO^-k^ 
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SBj 



ike 



20 



£2 
2 



SE'^, (52) 



In the ideal MHD limit, when SE'j^ = 0, one recovers 
magnetorotational mode. In the absence of rotation, a dis- 
persion relation for ideal MHD can be derived by setting 
determinant of left hand side matrix to zero. The departure 
from ideal MHD is due to the collisional effects. They will 
appear when electric field is eliminated in favour of magnetic 
field utilizing Ohm's law and Maxwell's equation. 

Making use of linearized Ampere's law 



SJ± = 



ike 
4tt 



1 
-1 



SB, 



in the generalized Ohm's law. 



SE', = 



-ike 



SOHr —up 



4 7rA(a;) \ up saui, 



SBj 



(53) 



(54) 



where A(a;) = aHr{io) aH<i>{ijj) + upiuj) and, s — sign{Bz). 
Introducing normalized variable y — i lu/Q, frequency depen- 
dent part of the conductivities can be written as. 



H.{y) = 
H2{y) = 



Q{y)^{l+PhF{y), 
(1 + Pl) Cl^ F{y) 



il+(3f)n2Fiy) 



(55) 



where F(jy) = Di{y)/D2{y) and Di{y) = 1 + ny,D2 = 
Df +Cli Cl2, y = iuj/n. 

Eliminating 6E'^ from equation (52), one gets the fol- 
lowing dispersion relation 



kvA 



(50) a{^\ +6 ^ +c = 



kVA 



(56) 



■ X G{yf + X G{y) (2 saH4> + 0.5 s ght + 2yop)a^^ 

+ (j/^ + 1) G{y), (57) 



b={2yGp-l.^sGH4,)(y +l)xG{y)(T^ 
+x'G(y)^(2j/^-3), 

c = X^G{yfy^{y^ + l). 



(58) 
(59) 



Here G{y) = l^{y)/a\, a±_ = \/op+^^ ^.nd the pa- 
rameter X = i^c/n is the normalized critical frequency 
ojc = 4 7r (va/c) a± similar to W99. As has been noted in 
W99, the ideal MHD description is valid in the large x limit. 
When kvA '^ ^c non- ideal MHD effect become important. 

When ion inertial effects are ignored, the expressions 
for a, bandc, in equations (57)-(59) reduces to W99. Ideal 
MHD is recovered in x — > co limit. As has been discussed 
in W99, the Hall term has considerable effect on the mag- 
netorotational instability growth rate (~ fJ) in the small x 
limit. This is because the Hall effect and collisions are in- 
tricately linked in a partially ionized plasma. Since the Hall 
effect is due to ion-neutral collision, small wavelength fluctu- 
ations are supressed in the Hall regime and only long wave- 
length fluctuations will grow. We shall see that both Hall 
as well as Ambipolar diffusion will become important in the 
small x, high frequency limit, suggesting that the presence 
of ion inertial effect destabilizes the weakly ionized disc at 
all wavelengths. 



5 ENERGETICS OF THE DISC 

Before we discuss numerical results, lets examine the various 
factors (viz. Lorentz force. Joule heating) that may affect 
magnetorotational instability . The Lorentz force, J X B, 
acts on the neutrals through collisions. Making use of equa- 
tion (31), it can be written as 



= (XpVb + XhB X Vi 



(60) 



J X B 

cpQ, 

Here xp = {<^p/<^±)x, Xh = {<^H/cri_)x and Vp = 
cE' X B /B is the drift velocity of B through the neu- 
trals. The first term on the right hand side is a measure 
of simultaneous acceleration and frictional drag; viewd from 
a neutral frame, this force accelerates the neutral towards 
E' X B velocity. The parameter xp provides the strength of 
the collisional coupling. With decreasing xp, i-e. when the 
ionized medium is far from ideal MHD regime, this term 
may become increasingly unimportant. Thus the modifica- 
tion to the ideal MHD modes will be severe in the small xp 
limit since collision modifies the fluid response to the mag- 
netic tension. Note that in the ambipolar regime, this term 
is responsible for dissipation as well as feeding of energy to 
the neutrals. The second term will accelerate the neutral in 
the direction of £?'. 

In order to understand the implications for the magne- 
torotational instability , we need to identify the conditions 
under which energy is fed to the fluctuations. Recall that 
the electric field E' = E + Vn X B/c is given in the neutral 
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frame and thus, 



J E = 



\\Eu + crp 



£;!' + «„•£! X B/( 



B 



— O-H [Vn ■ E'^) 



(61) 



Clearly then, the energy exchange consists of Joule heat- 
ing and acceleration of the neutral medium. The term 
o-||i5|l -\- (jpE'j_ is the Joule heating. This term is always 
positive for positive ctp . However, since ap may become neg- 
ative near resonance and the Ohmic term ap E'j_ may feed 
rather than dissipate energy. Therefore, in the ambipolar 
regime, fluctuations may grow. The terms w„ • E'j_ X B and 
Vn ■ E'j_, for ambipolar and Hall respectively, corresponds 
to the feeding or, extraction of the kinetic energy by the 
Lorentz force. Therefore, depending upon the sign of the ki- 
netic energy terms, the Lorentz force may facilitate either 
growth or damping of the magnetorotational instability . 



6 RESULTS 

We shall solve the dispersion relation (56) numerically in 
various limiting cases and discuss modifications due to ion 
inertia. In the absence of ion inertia, various x limits and 
its effect on the magnetorotational instability have been dis- 
cussed in detail in W99. We assume that electrons are frozen 
in the magnetic field, i.e. /3e — oo. In this limit, cr%/a% = Pi. 
Therefore, we shall solve the dispersion relation (56) by vary- 
ing key parameters Xi Pi and v = VinjO.. In accretion disc 
environment, the value of /3i, may vary in a wide range, we 
shall choose fii between 0.1 and 1. Although higher value of 
Pi can be chosen, the growth rate of magnetorotational in- 
stability will be very small. The parameter Uin/Q, is similarly 
varied in a wide range. 



6.1 Variation of Pi for x ~ 01, 1^ = 1 

In Fig. 2(a) we plot the growth rate by varying Pi while keep- 
ing V = 1 and X = 0.1 fixed. With the decreasing Pi, the 
parameter window of magnetorotational instability extends 
towards short wavelength and the growth rate exceeds ideal 
MHD limit for Pi — 0.1. With the decrease in Pi when 
<^ci < T^in ~ ^, the mode grows upto 0.92 fi for Pi = 0.1. 
This can be possibly attributed to the fact that owing to 
the faster coUisional (and Keplerian) time scales in compar- 
ison with the gyration time, rotational free energy becomes 
available to the fluctuations at the coUisional time scale. 
For further decrease in the value of Pi to 0.01, the maxi- 
mum growth rate decreases. This behaviour indicates that 
if Pi is increased beyond some critical value, the Hall diffu- 
sion is dominated by the ambipolar diffusion. The growth 
rate for Pi — 1 curve is small. This regime correspond to 
oJci ~ i^in ~ ^, i.e. the rate of ion gyration and collision with 
the neutrals is comparable with the rotational frequency. 
With the increase of Pi, the growth rate decreases and disap- 
pears altogether for very large Pi. Thus, Hall effect (caused 
by the ion inertial and coUisional effects), starts dominating 
the ambipolar diffusion and the mode starts growing. Fur- 
ther decrease in Pi and increase in Hall conductivity reduces 
the growth rate to 0.82^. The coUisional effect weakens with 
decreased Pi and thus, the parameter window operates in 
both small and large wavelength regimes. 



6.2 Variation of x and m for fixed Pi 

In Fig. 2(b) the growth rate is given for varying x a-nd fixed 
Pi = 0.1, !/ = 1. With the decreasing x (i-e. increasing coUi- 
sional coupling) , the wavenumber window of the instability 
shifts towards long wavelength, consistent with the fact that 
non-ideal MHD effects start playing an increasingly impor- 
tant role for smaller x- The result is similar to W99. This 
result is also in agreement with Blaes & Balbus (1994). We 
see that with the decreasing x the growth rate remains un- 
changed. Only change is in the wavenumber window that 
shifts towards the long wavelength consistent with the Hall 
dominated result of W99. 

In fig. 2(c) we plot the growth rate for X — 1 and 
Pi =0.1 for varying i^. The results are along the expected 
line. With the increase in collision, the growth rate decreases 
due to increased role of dissipation and the parameter win- 
dow of instability shifts towards long wavelength. When the 
ion-neutral collision rate is comparable (or smaller) to the 
rotational frequency, the MRI is unaffected by the colli- 
sion. This is because the rate of dissipative loss of the en- 
ergy is comparable or slower than the rotational time scale 
f^ ~ f2~^. Thus dissipation stops affecting the growth rate 
and it saturates around ~ 0.75. Further decrease of i/ does 
not affect the growth rate. 



6.3 General limit with Pi — 1 

In figure 3(a) we plot the growth rate for the positive ori- 
entation of the magnetic field with respect to the rotation 
axis, i.e. s = 1. The results are similar to the known re- 
sults of W99. However, there are some interesting differences 
towards small x limit. Whereas, the wavenumber window 
of magnetorotational instability shifts towards longer wave- 
lengths in small x limit, the growth rate of instability is 
not very sensitive. The ion inertia is able to provide the 
free energy to the fluctuations that can counterbalance the 
dissipative losses in small x regime where collision coupling 
between ion and neutral is very strong. Thus, with the de- 
crease in the value of x, small wavelength fluctuations are 
all suppressed leaving large wavelength modes to grow at 
~0.82!:!. 

In figure 3(b) we plot the growth rate against Uin/Q,. 
We see that when Vin ^ fi i.e., when collision time tc = v^^ 
is smaller than the rotational time tr = fi^^, and the free 
energy is dissipated by the collision. The growth rate of fluc- 
tuations decreases. We see from the plot that with increas- 
ing v the magnetorotational instability growth rate reduces 
significantly. However, the growth rate becomes insensitive 
beyond u — 10. This indicates that in the large u limit, 
ion inertial effect, namely Hall effect starts cancelling dis- 
sipation and thus, growth rate becomes insensitive to any 
further increase of v. In the opposite limit, i.e. when z^ ^ 1 
(or, tr < tc), the energy available to the fluctuation is at 
the rotational time scale and, infrequent, slow ion-neutral 
collision is unable to influence the growth of the instability. 
At 1/ = 0.01, the growth rate becomes maximum ~ O.Qil 
and any further decrease in v do not change the growth rate 
significantly. 
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Figure 2. The magnetorotational instability growth rate 
for Up — 1 for varying (3^ and gh with fixed v — 1 
and X = 0.1 is plotted against the wavenumber. The 
number against the curve is the value of ion Hall j3i for 
Oh = 1, 10 and 100. In figure 2(b) we hold v fixed and 
give growth rate for various v. In Fig. 2(c) we hold x 
fixed and vary v. 
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Figure 3. As for Fig. 2 but for ap — an = 1 



6.4 Weak ambipolar limit with Pi = 0.1, s = —1 

In figure 4(a) and 4(b) we plot the growth rate for Bz < 
by varying x E^nd v respectively. In figure 4(a) the growth 
rate is slightly smaller than for corresponding case in Fig. 
3(a) for s=l. For x ^ 0.1 the growth rate is insensitive to 
the changing value of x implying that if collision frequency is 
an order of magnitude smaller than the Keplerian frequency, 
the ambipolar effects are entirely compensated by the Hall 
and the growth rate remains constant. The curves are sim- 
ilar to W99 except that the mode exists for much smaller 
value of X than was the case in W99. Also, the growth rate 
is larger. The growth rate increases with increasing x and 
attains maximum value for x = 10. If x is increased further, 
there is no change in the growth rate. The rate at which ro- 
tational energy becomes available and dissipation operates 
become comparable. For x > 1, the ideal MHD limit is ap- 
proached and thus the effect of collision diminishes. Thus 
for X = oo the growth rate is maximum. 

In figure 4(b) the growth rate is not very sensitive to 
change in i/ except when it becomes large. For large v the 
growth rate reduces in comparison with the small f values 
and parameter window extends towards smaller wavelength. 
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Figure 4. As for Fig. 2(b), 2(c) but for s = -1. 
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Figure 5. As for Fig. 2 but with s = — 1. 



6.5 General limit for s = 1 and /3i = 1 

In Fig. 5(a) we plot the growth rate for s = — 1 for varying 
X- The growth rate is insensitive when x < 1- Small X is a 
measure of departure from ideal MHD and we note that Hall 
effects play important role for x < 1- Since for s — ~1 the 
sign of the wave helicity fl ■ SB is negative since Alfven wave 
is propagating in the negative direction (??). Thus, the in- 
crease in the non-ideal effect, manifested through Hall terms, 
does not have any bearing on the growth rate. With the 
increase of x the magnetorotational instability growth rate 
approaches ideal limit. 

In Fig. 5(b) the variation of growth rate with v is given. 
For !/ ^ 1 there is no change in the growth rate. For v > 1 
the instability is damped due to dissipation. The sign of the 
helicity ensures that non ideal effect do not feed the en- 
ergy to the fluctuations. Thus we see the decrease in the 
growth rate with increasing v. Therefore, we see that in the 
general case when the ambipolar conductivity may change, 
the growth rate of the instability becomes larger than the 
rotation frequency and the instability operates at long wave- 
lengths. Thus, the ion-inertial effect introduces entirely new 
feature to the dynamics of a weakly magnetized disc. Not 
only it changes the length scale over which the instability 
can operate but also, how fast it can operate. These feature 



makes ion inertial effect very important for the application 
to the protostellar discs. 



7 APPLICATIONS 

The modification to the magnetorotational instability by ion 
inertia may have wide ranging application in the astrophys- 
ical discs. Before discussing the application of the results, 
we shall note that the parameters x, ^ and /3i are not inde- 
pendent but are related by the following equation 






(62) 



where v = Vin/0. is the normalized collision frequency and 
X is the measure of non-ideal MHD effects. Therefore, the 
choice of x, Pi and u constrains the ratio pi/ pn and hence, 
the level of fractional ionization. 

At densities relevant to cloud cores and protostellar 
discs (densities Js lO^^cm""^), grains are the dominant 
charge carriers and their presence can significantly alter the 
dynamics of the disc. The ionization fraction is strongly af- 
fected by the abundance and size distribution of the grains 
through the recombination process on the grain surface. 
Near the midplane of the PPDs grains are the main charged 



Ion dynamics and the magnetorotational instability in weakly-ionized discs 11 



constituent (Wardle & Ng 1999). We shall assume a Keple- 
rian frequency for the minimum mass solar nebula (0.1 Mq), 
Q ~ 10~ s~ . Then, from equations (15) and (18) we write 



n 



^ 10" al: 



IQii cm- 



(63) 



Taking m„/mi ~ lO^^'* — 10^^^, we see that Ui„ is compa- 
rable to the dynamical frequency and ion-inertia becomes 
important. The collision of the energetic electron with neu- 
trals or magnetorotational instability induced turbulent con- 
vective homogenization of the entire disc (Inutsuka & Sano 
2005) may allow the magnetic field to couple to the disc 
matter near the mid plane. Therefore, the ion-inertia may 
modify the parameter window of instability near the mid- 
plane of the disc. 

In AGNs, for example NGC 4258, a thin disc of 0.2 pc 
diameter, bound by a central mass of --^ 2.1 x lO'^ ^01 is ro- 
tating with a velocity 900 kms"^ (GreenhiU et al. 1995). The 
observed emission emanates from an annulus of inner radius 
~ 0.13 pc and outer radius of 0.25 pc. Taking R = 0.1 pc, we 
get Q ~ 10^^" s^^. Thus the ratio Vin/0, ~ n„/cm~^. Taking 
ionization fraction Xe — 10^'' (Menou & Quataert 2001) at 
0.1 pc, we see from equation (29) that Rca = 10^^ n„/cm^^. 
The neutral density n„ ~ 10^ cm~^ and thus, both ion- 
neutral i^in as well as neutral-ion i^ni collision frequencies 
are very large in comparision with the rotation frequency. 
The charged grains are negligible in such a disk since from 
Pi/pn = 10~^, we have ni/n„ ~ 10"'^'* for micron-sized 
grains. Therefore, the charged grains are abesnt in such a 
disk and grain inertia will have no effect on the disc dynam- 
ics. 

Given the uncertain nature of the disc size, if we assume 
a disc of 100 pc with a temperature gradient towards the 
outer edge of the disc then the inertial effects in such a disc 
will be due to the charged grains near the core of the disc 
and due to the lighter ions near the surface region of the 
disc. For Q. ~ 10"* s"^ (at 100 pc), u^„/D, ~ 10"^ n„/cm"^, 
Xe ^ 1, we see that in the surface region of an disc, both ion 
and neutral fluid will be affected by the magnetorotational 
instability . Thus, ion-inertial efi'ect may be important in 
exciting MHD turbulence in the whole disc. 

Cataclysmic Variables (CV) are close binary systems 
with a white dwarf accreting material from a Roche-lobe 
and a companion low mass main sequence star. The typi- 
cal orbital frequency of the CVs vary between Q, ~ 10^'^ — 



lO-'^s-V Then i/,„/r2 



10- 



10 '^) X n„/cm 



For 



n„ ~ (10^ - lO'') cm"^, Ut„/n ~ 0(1). The temperature 
in CVs may vary in a wide range and disc can be mod- 
eled either as a weakly ionized plasma (Gammie & Menou 
1998) or a completely ionized plasma (Saxton et al. 2005). 
Clearly, both ion and neutral inertial effect operate on an 
equal footing in CVs. 

The circumnuclear disc at the galactic centre has a 
typical constant rotation speed of 110 Km s~^ ((Genzel & 
Townes 1987)) between 2 to 4 pc. The corresponding rota- 
tional frequency at 2-4 pc is, fi ~ 10"^^ - 10"" s"^ Then 
the ratio between ion-neutral coUison to the Keplerian rora- 
tion frequency is 



(10^ - 10^) X n„, 



(64) 



for given Vin value (equation (17)). Hence at 2 — 4 pc, the 
ion inertial response time, f2"^ is thousand times slower than 



the coUisional momentum exchange time, iy~^ . Therefore the 
ability of ion inertia to modify the magnetorotational insta- 
bility at 2 pc is unclear although Rca ~ 1 for ni/rin ~ 10""^. 
At 100 pc where fl drops by two orders, i/i„ becomes compa- 
rable to the rotational frequency and intertial effect on the 
magnetorotational instability may become important. 



8 CONCLUSIONS 

The paper examines the role of ion inertial effect on the mag- 
netorotational instability in a weakly ionized, thin, magne- 
tized Keplerian disc. The vertical stratification and radial 
and azimuthal variations were neglected - an approxima- 
tion valid for the wavelengths small compared to the disc 
scaleheight. The conductivity tensor becomes time depen- 
dent in the presence of ion inertial terms. This may result 
in the conductivity becoming negative near resonance. Fur- 
ther, radial and azimuthal component of Hall conductivity 
will be difi'erent. The following results were found. 

(i) The conductivities in a weakly ionized gas is in gen- 
eral complex in the presence of ion inertia and may become 
negative near the resonance points. The magnetorotational 
instability gets significantly modified in the presence of time 
dependent conductivities. 

(ii) In weak ambipolar regime, the presence of ion- 
inertial efi'ect substantially modifies the behaviour of the 
instability in the non-ideal (x < 1) limit. The maximum 
growth rate in the Hall dominated regime is ~ 0.92 (in the 
units of Q) and large wavelength fiuctuations can grow due 
to the inertial effect. 

For a fixed x, the growth rate is maximum (~ 0.8) for 
jy — 1 and starts to decrease with increasing v. Further, 
the parameter window shifts towards longer wavelength with 
increasing v. 

(iii) When both ambipolar and hall diffusion are com- 
parable, the maximum growth rate is ~ 0.85 and 0.95 for 
fixed u and x respectively. The increase in the ambipolar 
diffusion causes the fixed x case to have larger growth rate 
than fixed v case. 

To summarize, it is a common feature that large scale 
fluctuations exhibit the maximum growth rate of the magne- 
torotational instability when the ion inertial effects are in- 
cluded in the dynamics. This may have important impli- 
cation on the onset of turbulence in weakly ionized discs. 
For example, in PPDs, when grains are dominant ions, the 
grain inertial effect will significantly modify the magnetoro- 
tational instability growth rate and thus, will play an im- 
portant role in the onset of hydromagnetic turbulence. In 
AGNs also, grain will play important role. In CVs grain will 
provide the ion inertia away from the surface of the dwarf 
novae whereas, lighter ionized elements will provide the in- 
ertial effect close to the surface of the disc. Therefore, in CV 
discs, inertial effect may be important all across the diak 
and inertia modified magnetorotational instability may ef- 
fect the whole disc. In circumnuclear discs, the effect of ion 
inertia may be important far away from the centre of the 
disc. All in all, ion inertia seems to play an important role 
on the onset of magnetorotational instability . 

Present work inevstigates the role of ion inertia in the 
presence of an axial magnetic magnetic field. It will be in- 
teresting to investigate the role of ion inertia on the mag- 
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netorotational instability for a more general field geometry, 
particularly in the context of profile independent destabiliz- 
ing feature of ambipolar diffusion (Kunz & Balbus 2004). 
We shall leave this problem for future consideration. 
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